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The  two-dimensional  Smoluchowski  equation  is  employed  to  study  the  effect  of  elongational 
perturbations  on  nematic  liquid  crystal  polymers  under  a  weak  shear.  We  use  the  multiscale 
asymptotic  analysis  to  show  that  (1)  when  the  elongational  perturbation  is  small  relative  to  the  weak 
shear,  the  orientational  probability  density  function  (pdf)  tumbles  periodically  only  in  an 
intermediate  range  of  polymer  concentration;  outside  this  intermediate  range  (i.e.,  for  very  small  and 
very  large  polymer  concentration)  the  orientational  pdf  converges  to  a  steady  state  and  there  is  no 
tumbling.  (2)  When  the  elongational  perturbation  is  about  20%  of  the  shear  rate  or  larger,  the 
intermediate  range  of  tumbling  disappears  and  the  orientational  pdf  always  converges  to  a  steady 
state  regardless  of  the  polymer  concentration.  Our  theoretical  predictions  are  consistent  with  various 
earlier  results  based  on  the  Leslie-Ericksen  theory  [C.  V.  Chaubal  and  L.  G.  Leal,  J.  Non-Newtonian 
Fluid  Mech.  82,  22  (1999)]  or  analogous  3D  numerical  simulations  [M.  G.  Forest,  R.  Zhou,  and  Q. 
Wang,  Phys.  Rev.  Lett.  93,  088301  (2004);  M.  G.  Forest,  Q.  Wang,  R.  Zhou,  and  E.  Choate,  J. 
Non-Newtonian  Fluid  Mech.  118,  17  (2004)].  ©  2007  American  Institute  of  Physics. 
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I.  INTRODUCTION 

Nematic  liquid  crystal  polymers  (LCP)  are  technologi¬ 
cally  important  materials.*’^  For  example,  nematic  LCPs 
have  been  spun  into  high  strength  fibers  which  were  used  as 
airbags  that  cushioned  the  landing  of  NASAs  highly  success¬ 
ful  missions  to  Mars.^  The  dynamics  of  nematic  LCPs  has 
been  widely  modeled  by  the  kinetic  Doi-Hess  theory."^'^  The 
basic  idea  behind  the  Doi-Hess  theory  is  to  treat  the  LCP  as 
a  suspension  of  rigid  rodlike  nematogenic  molecules  and  de¬ 
scribe  the  ensemble  with  an  orientational  probability  density 
function  (pdf).  The  orientational  pdf  evolves  according  to  a 
nonlinear  Smoluchowski  (Fokker-Planck)  equation  where 
the  forces  acting  on  the  rodlike  molecules  include  hydrody¬ 
namic,  Brownian,  and  intermolecular  forces.  The  Smolu¬ 
chowski  equation  provides  a  rich  variety  of  mathematical 
problems  and  it  has  been  investigated  theoretically  for  pure 
nematic  equilibria,*"'^  extensional  flow-induced  equilibria,*^ 
equilibria  of  dipolar  ensembles,  *"*"'*  effect  of  high"^  and 
weak  shear,*^  and  effect  of  coplanar  magnetic  field. *^  Nu¬ 
merical  simulations  of  the  Smoluchowski  equations  and  the 
Doi  tensor  model  have  been  carried  out  extensively,  for  ex¬ 
ample,  see  Refs.  20-32. 

The  material  properties  of  the  LCPs  are  highly  affected 
by  the  processing  flows.  It  is  practically  important  to  under¬ 
stand  flow  effects.  Marrucci  and  Maffettone  gave  the  first 

33—35 

elegant  two-dimensional  analysis  on  the  effect  of  shear. 
Later  Larson^"  carried  out  three-dimensional  numerical 
simulations  using  spherical  harmonic  expansions.  Larson’s 
results  confirmed  the  theoretical  predictions  of  Marrucci  and 
Maffettone.  In  1999  Chaubal  and  Leal  applied  a  Lagrang- 


ian  numerical  method  called  the  smoothed  particle  hydrody¬ 
namics  technique  to  explore  the  effect  of  flow  perturbations 
on  LCPs  under  a  simple  shear.  They  considered  elongational 
perturbations  to  a  simple  shear  to  make  it  slightly  more  ex¬ 
tensional  than  a  pure  shear  and  found  that  slight  perturba¬ 
tions  to  the  flow  field  can  arrest  periodic  tumbling  behavior. 
This  phenomenon  has  been  observed  for  both  low  and  high 
shear  rates  even  though  the  mechanism  is  different.  In  Ref. 
20,  numerical  simulations  were  carried  out  in  the  case  of 
finite  Peclet  number  where  the  Peclet  number  denotes  the 
shear  rate  nondimensionalized  with  respect  to  the  rotational 
diffusivity.  In  the  asymptotic  limit  of  small  Peclet  number, 
the  Leslie-Ericksen  theory  was  used  to  explore  deviations 
from  simple  shear.  By  introducing  an  effective  aspect  ratio 
and  plotting  the  critical  value  of  perturbation  (a)  as  a  func¬ 
tion  of  the  polymer  concentration  ([/),  it  was  revealed  that 
for  £>'>0.04,  the  orientational  director  flow-aligns  for  all  {/; 
for  a  <0.04,  there  is  an  interval  I7n,ax)  within  which 

the  orientational  director  displays  periodic  tumbling  behav¬ 
ior;  outside  of  the  interval  (C^in,  the  director  fiow- 

aligns  again.  For  U  less  than  the  isotropic-nematic  transition 
4.49,  the  director  flow-aligns  for  all  a^Q. 

Recently  Forest’s  group  did  extensive  research  on  vari¬ 
ous  solutions  of  the  3D  Smoluchowski  equations  under  dif¬ 
ferent  kinds  of  conditions. In  particular,  in  Ref.  36 
they  formulated  a  monodomain  correspondence  principle  of 
kinetic  and  mesoscopic  theory  for  arbitrary  aspect  ratio  nem¬ 
atic  polymers  in  general  linear  planar  flow.  From  the  prin¬ 
ciple,  the  monodomain  response  of  nematic  polymers  for  all 
linear  flows  in  the  plane  of  shear  is  identical  to  the  mon- 
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odomain  response  of  another  nematic  polymer  liquid  with  a 
renormalized  molecular  aspect  ratio  parameter  in  pure  simple 
shear  flow.  Later  they  extended  the  study  in  Ref.  36  to  track 
the  boundaries  of  the  chaotic  region  of  kinetic  theory  due  to 
finite  molecule  aspect  ratio  and  the  addition  of  straining  flow 
in  the  plane  of  shear.  It  was  predicted  that  sheared  chaotic 
response  persists  up  to  a  threshold  straining  flow  strength 
and  minimum  aspect  ratio,  beyond  which  chaotic  behavior  is 
arrested.  Both^^’^®  gave  predictions  on  how  steady  and  un¬ 
steady  attractors  are  modified  and  eventually  lost  or  persist  in 
the  presence  of  a  variable  extensional  flow  contribution. 

In  the  2D  case,  it  is  well-known  that  in  the  absence  of 
flow  the  isotropic-nematic  phase  transition  occurs  at  U=2, 
where  U  is  the  normalized  polymer  concentration.  In  the 
presence  of  an  imposed  weak  shear  there  is  a  threshold  {Uq 
—  2.41144646)  for  U:  When  U<Uq,  steady  state  solution 
exists;  otherwise  there  is  no  steady  state  and  the  orientational 
pdf  is  temporally  periodic  (“tumbling”)  and  can  never  reach 
a  steady  state.  If  one  perturbs  the  weak  shear  by  a  small 
elongation,  this  scenario  will  change.  In  fact,  a  shear  flow 
can  be  viewed  as  the  superposition  of  a  rigid  body  rotation 
and  an  elongational  flow.  So  the  imposed  elongational  flow 
will  either  increase  or  decrease  the  elongational  component 
of  the  shear  flow  depending  on  how  their  directions  are 
aligned  with  respect  to  each  other.  We  will  show  that  if  the 
elongational  component  is  increased  by  a  small  fraction,  then 
there  is  no  tumbling  at  large  polymer  concentration.  This 
leads  to  the  situation  that  the  orientational  pdf  tumbles  peri¬ 
odically  only  in  an  intermediate  range  of  polymer  concentra¬ 
tions;  outside  this  intermediate  range  (below  the  lower  end  of 
the  range  or  above  the  upper  end  of  the  range)  the  orienta¬ 
tional  pdf  converges  to  a  steady  state  and  there  is  no  tum¬ 
bling.  If  the  elongational  component  is  increased  further  (by 
about  20%  of  the  original  elongational  component  of  the 
shear  flow),  the  intermediate  range  of  tumbling  disappears; 
that  is,  for  all  polymer  concentrations,  the  orientational  pdf 
converges  to  a  steady  state.  Our  theoretical  results  are  con¬ 
sistent  with  Refs.  20,  25,  and  36. 

■3  0 

Our  asymptotic  analysis  invokes  a  multiscale  analysis 
with  the  assumption  that  the  solution  is  a  quasisteady  state. 
From  the  physical  point  of  view,  this  assumption  is  reason¬ 
able.  When  the  flow  field  is  weak,  it  takes  a  long  time  for  the 
flow  field  to  change  the  orientation  distribution  significantly. 
In  other  words,  the  time  scale  of  the  flow  field  affecting  the 
orientation  distribution  is  large.  On  the  other  hand,  in  the 
absence  of  flow  field,  the  time  scale  of  the  intermolecular 
force  (Maier-Saupe  interaction)  relaxing  the  orientation  dis¬ 
tribution  to  the  Boltzmann  form  is  fixed.  So  in  the  presence 
of  a  weak  flow  field,  it  is  reasonable  to  expect  that  the  inter¬ 
molecular  interaction  will  be  able  to  drive  the  orientation 
distribution  approximately  to  the  Boltzmann  form  before  the 
flow  field  can  take  a  significant  effect  on  the  orientation  dis¬ 
tribution.  This  is  because  the  time  scale  of  relaxing  to  the 
Boltzmann  form  is  small  in  comparison  to  the  time  scale  of 
the  weak  flow  field.  Thus,  it  is  reasonable  to  assume  math¬ 
ematically  that  to  the  leading  order  the  orientation  distribu¬ 
tion  is  of  the  Boltzmann  form.  Of  course,  the  director  of  the 
orientation  distribution  may  change  in  time  driven  by  the 
weak  flow  field.  Although  this  assumption  is  reasonable  on 


physical  grounds,  it  is  desirable  to  have  a  rigorous  math- 
ematical  justification  based  on  the  solvability  condition  or 
other  arguments.  Currently  such  a  rigorous  mathematical  jus¬ 
tification  remains  an  open  problem.  We  will  pursue  this  issue 
in  our  future  work.  In  recent  years,  multiscale  asymptotic 
analysis  has  become  a  powerful  mathematical  tool  in  many 
applications.  In  Ref.  39,  Choate  and  Forest  used  it  to  inves¬ 
tigate  the  viscoelastic  response  of  nematic  polymers  to  small 
amplitude  oscillatory  shear  using  the  Doi-Hess  mesoscopic 
orientation  tensor  model;  Vicente  Alonso  et  al."^  applied  it  to 
study  the  nonlinear  dynamics  of  a  nematic  liquid  crystal  in 
the  presence  of  a  plane  Couette  flow  using  a  Landau-de 
Gennes  model;  Chillingworth  et  employed  it  to  describe 
the  geometry  and  dynamics  of  a  nematic  liquid  crystal  in  a 
uniform  shear  flow. 

It  is  worthwhile  to  point  out  that  other  perturbation 
schemes  have  been  used  to  study  both  the  kinetic  and  tensor 
models  for  nematic  liquid  crystal  polymers.  For  example, 
Kuzuu  and  Doi^^’^^  and  Semenov"^^  developed  a  weak-flow 
asymptotic  analysis  of  the  molecular  kinetic  equation  for 
nematic  liquid  crystals,  which  provides  a  molecular  basis  for 
all  continuum  theory  parameters.  Doi  and  Larson^^  analyzed 
the  effect  of  steady  flow  fields  on  the  isotropic-nematic  phase 
transition  of  rigid  rodlike  polymers  by  a  singular  perturba- 
tion  method;  Forest,  Zhou,  and  Wang  expanded  the  orien¬ 
tational  distribution  function  in  the  Peclet  number  to  find 
approximate  steady  solutions  in  weak  shear.  The  works  of 
Refs.  42,  43,  45,  and  27  all  solved  the  kinetic  equation  and 
derived  rigorous  “solvability  conditions.”  The  advantage  of 
this  kind  of  approach  is  that  unlike  multiscale  asymptotic 
analysis,  it  does  not  posit  a  form  of  the  solution.  However, 
these  earlier  works  focused  on  the  persistence  of  steady 
states  and  did  not  perform  constructions  of  asymptotic  peri¬ 
odic  solutions.  Previous  results  on  periodic  solutions  have 
been  given  by  Lee,  Forest,  and  Zhou^®  using  a  mesoscopic 
2D  tensor  model.  By  invoking  the  Poincare-Bendixson  theo¬ 
rem  which  describes  the  long  term  behavior  of  orbits  of  con¬ 
tinuous  dynamical  systems  on  the  plane,  they  gave  an  el¬ 
egantly  simple  proof  that  limit  cycles  (“tumbling  orbits”) 
must  arise  beyond  the  parameter  boundary  for  the  steady- 
unsteady  transition. 

The  purpose  of  this  paper  is  to  revisit  the  well-studied 
elongation-perturbed  shear  problem  from  a  different  theoret¬ 
ical  point  of  view.  Our  study  will  provide  a  detailed  analysis 
on  how  perturbations  of  flow  conditions  affect  the  system 
dynamics  at  small  flow  rates.  We  will  develop  a  multiscale 
analysis  directly  on  the  Smoluchowski  equation  of  the  ki¬ 
netic  theory.  We  thus  extend  results  of  Chillingworth  et 
Vicente  Alfonso  et  Lee  et  and  Choate  et  from 
a  tensor  model  to  the  kinetic  model.  Our  work  differs  from 
Refs.  25  and  36  in  three  parts:  (1)  Refs.  25  and  36  are  based 
on  numerical  simulations  of  the  3D  Smoluchowski  equation, 
whereas  ours  is  on  the  asymptotic  analysis  of  the  2D  Smolu¬ 
chowski  equation;  (2)  Refs.  25  and  36  include  the  effect  of 
arbitrary  aspect  ratio,  whereas  our  study  is  limited  to  the 
infinite  aspect  ratio;  and  (3)  Refs.  25  and  36  concentrated  on 
the  persistence  of  steady  states,  whereas  our  focus  is  on  both 
steady  and  unsteady  solutions. 

In  this  paper  we  limit  our  attention  to  the  two- 
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dimensional  case  with  a  modified  weak  shear  flow.  More 
precisely,  the  Peclet  number  from  shear  flow  (e)  is  small  and 
the  Peclet  number  from  elongation  {q  ■  e)  is  of  the  same  order 
or  smaller  than  that  of  the  shear  flow.  Our  restriction  of  the 
orientational  probability  density  from  a  function  on  the 
sphere  to  the  circle  is  mainly  due  to  its  mathematical  sim¬ 
plicity,  its  possible  relevance  to  experimental  results  on 
monolayer  films  of  nematic  polymers,"^^’"^^  and  its  potential 
guidance  to  the  full  3D  orientational  distributions.  This  kind 
of  approach  has  been  widely  adopted  by  others. 

This  paper  is  organized  as  follows:  We  give  a  brief  de¬ 
scription  on  the  mathematical  formulation  of  the  Smolu- 
chowski  equation  coupled  with  modified  shear  flow  for  nem¬ 
atic  liquid  crystalline  polymers  in  Sec.  II.  We  show  the  effect 
of  elongational  perturbations  on  the  weak  shear  behavior  in 
Sec.  III.  Finally,  in  Sec.  IV  we  summarize  our  results. 


II.  MATHEMATICAL  FORMULATION 
FOR  THE  TWO-DIMENSIONAL  CASE 


y  y 

v=-(y,-x)-F-(-y,-x)  =  v„H-V£, 


(5) 


where  y  is  the  shear  rate.  For  the  elongational  flow 
=  y/2(-y  ,-jr),  we  have 


-1  0  J 
V2  V2 

V2  V2 

2  2 

1  0 

2  ~  2 

V2  V2 

.0  -  1  . 

V2  V2 

~  2  2 

2  2 

(6) 


So  Vvg  has  eigenvector  [V2/2,-v5/2]^  associated  with  the 
eigenvalue  1  and  eigenvector  [^272,^2/2]^  associated  with 
the  eigenvalue  -1. 

The  torque  induced  by  the  rotational  flow  on  a  poly¬ 
mer  rod  is 


We  briefly  review  the  two-dimensional  mathematical 
formulation  of  the  Doi-Hess  kinetic  theory  for  homogeneous 
flows  of  rigid  rodlike  nematogenic  molecules  immersed  in  a 
viscous  solvent  subject  to  an  imposed  flow  field."*’^  First,  we 
denote  the  orientation  of  each  polymer  rod  by  an  angle  9. 
Then  the  orientational  direction  of  a  polymer  rod  is  given  by 
u=(cos  0,sin  0).  To  model  the  nematic  tendency  of  the  rod¬ 
like  molecules  to  align  with  each  other,  it  is  common  to  use 
the  short  range  Maier-Saupe  potential, 

Vms  =  “  2A:5ri7uu:(uu) 


kgT  y  e 
.^^  =  ksT-, 
Dr  2  "2 


(7) 


where  Dr  is  the  rotational  diffusivity  of  the  polymer  rod 
(assumed  constant  here),  and  e=  yIDr  is  the  Peclet  number 
of  the  shear  flow.  Similarly,  the  torque  induced  by  the  elon¬ 
gational  flow  V£  is 


kgT  y  B 

- cos  29=  -  koT—  cos  29. 

Dr  2  "2 


(8) 


=  -  kgTU[{cos  20)cos  29+  (sin  29)sm  29], 
where  uu  and  A :  B  are  defined  as 


(1)  If  we  increase  the  elongation  by  a  fraction  of  q,  then  the  total 
elongation  is 


/  UiUi  U1U2  \ 

uu  ^ 

\M2M1  M2M2/ 

(2) 

A:B  =  2  cupij. 
ij 

In  Eq.  (1)  17  is  the  normalized  polymer  concentration  which 
describes  the  strength  of  intermolecular  interactions,  kg  is  the 
Boltzmann  constant,  T  is  the  absolute  temperature,  and  (uu) 
is  the  second  moment  of  the  orientation  distribution, 

(uu)  =  I  uup(u,x,t)<7u,  (3) 

||u||=l 

where  p(u,x,f)  is  the  orientational  probability  density  func¬ 
tion  of  the  ensemble,  i.e.,  the  probability  density  that  a  poly¬ 
mer  rod  has  direction  u  at  location  x  and  time  t. 

A  shear  flow  with  velocity 

v=y(0,-x)  (4) 


(9) 

and  the  torque  induced  by  the  total  elongation  is 

F£j.=  -^b7’“(1 +?)cos26>.  (10) 

Figure  1  depicts  velocity  fields  of  a  shear  flow  and  an 
elongation-perturbed  shear  flow,  respectively. 

If  we  assume  the  polymer  orientation  distribution  is  ho¬ 
mogeneous  in  space,  then  the  orientational  pdf  under  an  im¬ 
posed  perturbed  shear  evolves  according  to  the  Smolo- 
chowski  equation,'^ 

where  Vms(^)  the  Maier-Saupe  potential  (1)  and  yps(^)  is 
the  potential  due  to  the  external  field  (i.e.,  perturbed  shear 
flow) 


can  be  decomposed  into  the  sum  of  a  rotational  flow  and  y  )  sin  20 

an  elongational  flow  Vg,  ^^24^ 


(12) 
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Here,  for  simplicity,  all  potentials  have  been  normalized  by 
kgT.  In  Eq.  (11),  we  can  set  D^=  1  by  introducing  the  nondi- 
mensionalized  time  Below  we  are  going  to  do  mul¬ 

tiscale  asymptotic  analysis  on  (11)  with  £>,.=  1. 


III.  MULTISCALE  ASYMPTOTIC  ANALYSIS 


We  consider  Eq.  (11)  for  small  e  which  allows  us  to 
carry  out  a  multiscale  asymptotic  analysis.  First,  we  simplify 
the  expression  of  the  Maier-Saupe  potential  (1)  by  selecting 
an  angle  a  so  that  (sin  2(0-0'))  =  0.  Then  the  Maier-Saupe 
potential  can  be  written  in  a  compact  form, 

^Ms(^)  =  “  f/(cos  2(0-  Q!))cos  2(0-  a).  (13) 


We  restrict  our  application  of  the  multiple  scale  method 
to  the  introduction  of  two  time  scales,  Tf^  =  t  and  Ti  =  e-t. 
Assume  that  p  possesses  an  expansion  of  the  form 

p(0,f)  =  p^°\e,To,Ti)  +  8p(')(0,ro,ri)  +  oi^^)  (i4) 


and  substitute  it  into  Eq.  (11)  to  obtain  the  leading  order 
equation 


dp' 


,(o) 


07'n  dO 


(15) 


Here  V^j(0)  =  -[/(cos  2(0-0))*°'  cos  2(0-a)  corresponds  to 
the  leading  order  term  of  the  Maier-Saupe  potential  and  the 
average  in  the  expression  is  taken  with  respect  to  the  pdf 


p*°'(0).  Since  Eq.  (15)  does  not  contain  an  external  field  ef¬ 
fect,  we  have  a  pure  nematic  equation  at  the  leading  order. 
Also  notice  that  Eq.  (15)  does  not  contain  time  scale  Tj.  As  a 
result,  the  orientation  distribution  will  converge  to  a  steady 
state  with  respect  to  time  scale  Tq  (it  may  not  be  a  steady 
state  with  respect  to  time  scale  Tj).  The  steady  state  pdf  of 
Eq.  (15)  is  given  by  the  Boltzmann  distribution. 


p*°'(0,7’o,ri)  = 


exp[r  cos  2(0-  a)] 
/o'^exp[r  cos  2(0-  a')](f0’ 


(16) 


where  r  is  defined  as 

r=  U{cos  2(0-a))*°'. 


(17) 


Equation  (17)  is  also  the  equation  for  determining  r  for  a 
given  value  of  normalized  polymer  concentration  U  in  the 
absence  of  flow.  In  the  case  of  weak  shear  plus  weak  elon- 
gational  perturbation,  the  leading  term  of  r  satisfies  Eq.  (17). 
We  consider  the  simple  situation  where  the  polymer  orienta¬ 
tion  starts  with  the  steady  state  pdf  (16).  Then  Eq.  (15)  im¬ 
plies  that  the  pdf  remains  unchanged  with  respect  to  Tq. 
Mathematically,  this  means 


0p*°'(0,7’o,ri) 


0Tn 


=  0. 


(18) 


So  p^°\e,To,Ti)  is  independent  of  Tq  and  thereby  “quasi¬ 
steady.”  As  a  result,  the  phase  angle  a  is  also  independent  of 
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Tg-  But  a  may  depend  on  and  this  corresponds  to  the  slow 
time  evolution  invoked  by  the  imposed  flow  field.  To  find 
q'(7’i),  consider  the  expansion  of  dp{0,t)ldt, 


dp{e,i)  dp^^Kejoj,) 


dt 


dT. 


+  e 


dp^°\e,To,T,) 


dT, 


dp^\0,TQ,T{) 


dTa 


(19) 


Both  our  extensive  numerical  simulations  and  the  physical 
intuitions  suggest  that  p^^\B,Tq,T\)  is  independent  of  time 
scale  Tq.  Although  a  rigorous  mathematical  proof  is  still  an 
open  problem,  here  in  our  multiscale  expansion,  we  adopt 
the  assumption  that  p^^\6,Tq,T{)  is  quasisteady.  With  this 
assumption,  Eq.  (19)  becomes 


dp{e,t)  dp(^\e,TQ,T{) 


+  0{e^) 


dt  dTi 

JP^^\0,Tq,T,)  da{T,) 
^  00  '  dT, 


+  0{e^). 


(20) 


On  substituting  this  into  Eq.  (11)  and  rearranging  terms  we 
obtain 


0 

00 


^  +  ^(1  +  ^)cos  20+  VU0)jp{e)  +  ^ 


H- -H  0(8^) 


dT, 


=  0. 


(21) 


This  suggests 


■  ^  +  ^(1  +  ^f)cos  20+  yMs(^) 


p(p)  +  M^  +  ep(0) 


00 


X(0)^^  +  O(82)  =  -e/i, 

ail 


(22) 


where  eJi  is  the  integration  constant  of  the  order  (9(e).  In 
Eq.  (22)  both  da{T-i)ldTi  and  are  unknown  so  we  need 
two  linear  equations  to  determine  them.  To  do  so,  we  first 
integrate  Eq.  (22)  from  0  to  277  to  obtain 

fllT  rlTT 

V^^{0)p{0)d0 

Jo 


---H-(lH-g)^  cos20p{0)d0  + 


2  2 


da{Ti)  , 

+  8 -  +  0{e^)  =  —  2778/[  . 

dT^ 


(23) 


The  second  term  on  the  left  of  Eq.  (23)  has  the  expression 

r2'7r 

I  cos20p{0)d0={cos20) 

Jo 

=  (cos  2(0-  a)cos  2a 
-  sin  2(0-  Q;)sin  2a) 

=  (cos  2(0-  a)cos  2a) 

=  (cos  2(0-  a))^°^  cos  2a  +  0(e) 


=  —  cos  2a  +  0(e), 
U 


(24) 


where  rlU=  (cos  2(  0-  is  the  order  parameter.  The 

third  term  on  the  left  of  Eq.  (23)  has  the  expression 

rliT  rliT  ^ 

v;^si(^)p(0)d0=  — l-r  cos  2(0- a)]p(0)d0 

Jo  Jo 


-f 

Jo 


2r  sin  2(0-  a)p(0)d0 


■■  2r(sin  2(0-  a))  =  0, 


(25) 


where  we  have  used  the  fact  that  the  angle  a  is  selected  to 
make  (sin  2(0- a)) =0.  With  these  expressions,  Eq.  (23)  be¬ 
comes 

1  .  i+q  r  ,  da(Ti)  ^ 

-  -  +  COs[2Q'(ri)]  -F  +  0(8)  =  -  277/1. 

2  2  U  dll 

(26) 

To  find  a  second  equation  for  da(Ti)ldTi  and  Ji,  we  ma¬ 
nipulate  Eq.  (22)  by  dividing  it  by  ep(0),  approximating  p(0) 
by  p^^\0),  and  integrating  from  0  to  277  to  obtain 

[  (--)d0+  [  -(l+q)cos20d0+- f  V^sdB 

Jo  \  2/  Jo  2  eJo 

1  0p<°^(0)  p^da(Ti) 


p(°)(0)  00 


-d0  + 


f 

Jo 


dT, 


-d0+  0(e) 


P^°\0) 


d0. 


(27) 


Note  that  the  integrands  cos  20,  and  H p^^\0)0p^^^ 

X(0)/00=0 ln[p^°*(0)]/00  are  all  periodic  on  [0,27r],  so  the 
second,  third,  and  fourth  definite  integrals  in  Eq.  (27)  vanish. 
Using  Eq.  (16)  and  carrying  out  the  details,  we  get 


-  '^+2’Tr^  j  H-  0(e)=-(2TT)^g(r)Ji, 
dl  I 

where  the  function  g(r)  is  defined  as 


§('•)  = 


f 

Jo 


—  I  exp(r  cos  20)c/0 
277 


(28) 


(29) 


Multiplying  Eq.  (26)  by  a  factor  27rg(r)  and  subtracting 
it  from  Eq.  (28)  leads  to 


da(Ti) 

77-F277— —  =  277g(r) 

di  1 


-  + - COS  2a(ri) 

2  2  U 


da(T]) 

+ 

dTy  j 

Solving  for  da(T {} ! dTi  yields 


(30) 


da(Ty)  1 

dT,  “2[g(r)-l] 

Xcos[2a'(ri)] 


1  H-g(r)  1 


(l+q)r 

U 


(31) 


Applying  the  double  angle  formula  cos[2Q'(ri)]=  1 
-2  sin^[a'(7’i)],  we  obtain  after  some  manipulations 
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da{Ti) 

dTi 


=  Cl 


•{sin2[Q!(7’i)]  +  C2}, 


(32) 


where  Ci  and  C2  are 

gir) 


Cl- 


Cl'- 


gir)  -  1  U 
-  1 


“  •  (!+?)> 


(33) 


-[pir)  +  q]. 


2(1 +q)' 

Here  the  function  pir)  is  defined  as 


Pir)^-  ^ 
r  L  gin 


1 


H-  1. 


(34) 


Note  that  in  the  order  parameter  rlU={cos2i9-a)}^^\ 
probability  density  (16)  is  used  in  calculating  the  average.  As 
a  function  of  r,  r!  U  has  the  expression 

r  cos  20  exp(r  cos  261)^0 

U  /Q'^exp(r  cos  20)cf0 


(35) 


As  a  matter  of  fact  Eq.  (35)  is  the  equation  for  determining  r 
for  a  given  value  of  normalized  polymer  concentration  U. 
Here  for  mathematical  convenience,  we  use  r  (instead  of  U) 
as  the  independent  variable.  Note  that,  for  r>0,  the  mapping 
between  r  and  U  as  given  in  Eq.  (35)  is  monotonic  (see 
Appendix  A  for  a  mathematical  derivation  of  the  monotonic¬ 
ity).  As  a  consequence  of  the  monotonicity,  for  r>0,  the 
mapping  between  r  and  U  is  one-to-one.  Therefore,  a  result 
in  terms  of  r  can  easily  be  translated  to  the  corresponding 
result  in  terms  of  U.  Below  we  will  first  discuss  the  dynamic 
behavior  of  q;(7’i)  governed  by  Eq.  (32)  for  various  regions 
of  r.  Then  all  conclusions  will  be  restated  in  terms  of  U,  the 
normalized  polymer  concentration. 

Equation  (32)  describes  the  dynamics  of  the  slow  time 
evolution  of  the  phase  angle  q;(7’i).  In  Eq.  (32),  the  factor  ci 
is  always  positive  for  r>0.  If  C2  is  negative  (and  |c2|<l), 
then  the  phase  angle  aiT\)  will  stop  at  an  angle  ag  satisfying 
sin^  a'oH-C2=0.  In  other  words,  the  polymer  orientation  dis¬ 
tribution  will  converge  to  a  steady  state.  If  C2  is  positive,  then 
the  phase  angle  aiT\)  will  keep  increasing  and  the  polymer 
orientation  distribution  will  rotate  periodically  in  time  (tum¬ 
bling).  Thus,  positive  C2  causes  tumbling.  The  sign  of  C2  is 
the  opposite  of  the  sign  of  pir)  +  q.  It  follows  that  negative 
value  of  pir)  +  q  corresponds  to  tumbling.  Eigure  2  shows  the 
graph  of  function  pir)  and  its  asymptotic  approximation  as 

r^  H-oo, 

e('-)  =  ^(i  +  ^)+  as  H-oo.  (36) 


The  asymptotic  approximation  of  pir)  as  H-oo  is  derived 
in  Appendix  B.  As  shown  in  Fig.  2,  the  function  pir)  starts  at 
r=0  with  p(0)  =  1.  As  r  increases  from  0  to  r^^^,  the  function 
pir)  decreases  monotonically.  In  particular,  the  function  pir) 
intersects  with  the  x-axis  at  ro=  1.3276.  At  rnjin=2.6463,  the 
function  pir)  attains  its  minimum,  and  the  minimum  value  is 
min^p(r)  = -0.1 8754.  For  r>rinin,  function  pir)  increases 
monotonically.  For  r>rjni„,  function  pir)  is  well  approxi- 


FIG.  2.  (Color  online)  Graph  of  function  p(r)  and  its  asymptotic  approxi¬ 
mation  as  r— >  -l-“.  The  function  p{r)  intersects  with  the  x-axis  at  Tq 
=  1.3276.  At  r„,i„=2.6463,  the  function  p{r)  attains  its  minimum,  and  the 
minimum  value  is  min,.p(r)=-0. 18754. 


mated  by  its  asymptotic  expression  at  r=+^. 

Figure  3(a)  shows  the  tumbling  region  (shaded)  in  the 
iq,r)  space.  Here  the  tumbling  region  is  the  largest  region 
where  pir)+q  is  negative.  Recall  from  our  discussion  above 
that  negative  value  of  pir)  +  q  corresponds  to  tumbling.  So  if 
the  value  of  iq ,  r)  falls  in  the  tumbling  region,  then  the  poly¬ 
mer  orientation  rotates  periodically  in  time  (tumbling).  From 
the  asymptotic  expression  of  function  pir)  for  large  r,  we  see 
that  for  small  q,  the  upper  boundary  of  the  tumbling  region 
in  the  iq,r)  space  is  approximated  by 


2r 


1  /  3 

0'  =  :r  1  +  7“ 


4r 


(37) 


Solving  r  in  terms  of  q,  we  obtain 
1  H-  Vl  -I-  6q 

'•upper  =  —  ■ 


(38) 


This  asymptotic  expression  is  valid  for  small  q. 

Because  of  the  one-to-one  correspondence  between  r 
and  U,  a  region  in  the  iq ,  r)  space  can  be  uniquely  mapped  to 
a  region  in  the  iq,U)  space.  Figure  3(b)  shows  the  corre¬ 
sponding  tumbling  region  in  the  iq ,  U)  space,  where  U  is  the 
normalized  polymer  concentration  and  q  is  the  magnitude  of 
the  elongational  perturbation  relative  to  the  original  elonga- 
tional  component  of  the  shear  flow.  Again,  for  small  q  (cor¬ 
responding  to  large  r  and  large  U),  the  upper  boundary  of  the 
tumbling  region  in  the  iq ,  U)  space  can  be  well  approxi¬ 
mated  by  using  the  asymptotic  expression  of  function  piU) 
for  large  U  (which  is  derived  and  presented  in  Appendix  B). 
For  small  q,  the  upper  boundary  in  U  satisfies  approximately 


^  2U\  4Uj 


(39) 


Solving  U  in  terms  of  q,  we  obtain 
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FIG.  3.  (Color  online)  Tumbling  regions  are  shown  as  shaded,  (a)  Tumbling  region  in  the  {q,r)  space,  (b)  Tumbling  region  in  the  {q,U)  space.  The  critical 
value  for  q  is  ^^,=0. 18754. 


1  H-  Vl  H-  lOfl 

u  = - - 

Supper  4^ 


(40) 


This  asymptotic  expression  is  valid  for  small  q.  In  Fig.  3,  the 
asymptotic  approximation  of  the  upper  bound  of  the  tum¬ 
bling  region  for  small  q  is  shown  as  a  dashed  line  and  it 
agrees  well  with  the  exact  solution. 

Let  g^  =  -min^/7(r)=0. 18154.  q^  is  the  critical  value  of 
q,  when  the  q>qc  function  p{r)+q  is  positive  for  all  values 
of  r>0;  when  0<q<qc  the  function  p{r)+q  may  be  nega¬ 
tive  for  certain  values  of  r.  As  shown  in  Fig.  3,  when  0<q 
<q^  tumbling  occurs  only  for  r  (U)  in  a  finite  interval.  As  q 
increases  toward  the  critical  value  q^,  this  tumbling  interval 
gets  smaller  and  smaller.  When  q  >  q^,  the  tumbling  interval 
disappears  and  there  is  no  tumbling.  That  is,  with  an  elonga¬ 
tional  perturbation  of  less  than  20%  (^^=0. 18154)  of  the 
original  elongational  component  of  the  shear  flow,  the  tum¬ 
bling  disappears  for  all  values  of  U. 
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APPENDIX  A:  MONOTONICITY  OF  U(r) 

We  first  apply  integration  by  parts  to  the  numerator  on 
the  right-hand  side  of  Eq.  (35), 


cos  2^exp(r  cos  26)d0 


IV.  CONCLUSIONS 

We  have  applied  the  multiscale  asymptotic  analysis  to 
study  the  effect  of  elongational  perturbations  on  the  weak 
shear  behavior  of  nematic  liquid  crystal  polymers.  It  is  found 
that  when  the  elongational  perturbation  is  small,  the  tum¬ 
bling  behavior  occurs  only  in  an  intermediate  range  of  poly¬ 
mer  concentrations;  when  the  elongational  perturbation  is 
above  about  20%  of  the  shear  rate,  tumbling  disappears  for 
all  polymer  concentrations.  These  findings  suggest  a  mecha¬ 
nism  on  the  control  of  the  structure  of  nematic  liquid  crystal 
polymers. 


=  -  I  exp(r  cos  26)d{sm  20) 

2  Jo 

rlTT 

=  r\  sm^20exp{rcos20)d0 

Jo 

rliT 

=  r\  (1 -cos^  20)exp(rcos  20)(i0.  (Al) 

Jo 

Substituting  this  result  into  Eq.  (35)  and  writing  {/  as  a  func¬ 
tion  of  r,  we  have 
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1  /q’^cos^  20exp(r  cos 


U{r) 


/o’^exp(r  cos  9)d9 


(A2) 


f 

Jo 


exp(r  cos  29)d9 


To  show  that  U{r)  is  monotonically  increasing,  we  only  need 
to  show 


=/' 


■■  2  I  exp[r(l  -  2s)]j=-i^=ds 
>0  isil-s 


fo^cos^  29exp(rcos29)d9 
/o’^exp(r  cos  29)d9 


(A3) 


Jo 


:  2  exp(r)  I  exp(- 2^^)^|  1 -I- -I- 


ds 


is  monotonically  increasing.  Function  H'(r)  can  be  viewed  as 
an  average,  H'(r)=(cos^  29)  with  probability  density 


p{9,r): 


exp(r  cos  2  9) 

/o'^ exp(r  cos  29)d9 


(A4) 


The  derivative  of  the  probability  density  with  respect  to  r  is 
dp{9,r) 


^  ,Jr(i/2)  1  r(3/2) 


,,,  11 
^expW^l  +  -.-  + 


Substituting  this  result  into  function  g{r),  we  have 


(B3) 


dr 


:  (cos  2  0  -  (cos  2  9))  pi  9,  r) . 


(A5) 


■iU 


gir)  =  —  I  exp(r  cos  29)d9 


27r  r 


(A6) 


That  is,  as  -Foo,  function  gir)  is  exponentially  large  and 
consequently  1  /gir)  is  exponentially  small.  In  a  similar  way, 
we  obtain 


rlTT 

J„  ' 


(1  -  cos  26)exp(r  cos  29)d9 


ds 


The  hrst  derivative  of  w(r)  is 

w'ir)  =  (cos^  2^(cos  29- {cos  29))). 

The  second  derivative  of  w(r)  is 

w"ir)  =  (cos^  20(cos  2^-  (cos  29))^) 

-  (cos^  29)  ■  ((cos^  29)-  (cos  29)^) 

=  - 2(cos  20  w'(r) -I- var(cos^  26) .  (A7) 

Thus,  function  wir)  satishes  the  special  property  that  when¬ 
ever  w'(r)=0  we  have  w"ir)>Q.  At  r=0,  we  have  >^(0) 

=  0.5  and  w'(0)  =  0.  The  special  property  of  wir)  guarantees 
that  w'(r)>0  for  r>0.  This  leads  immediately  to  the  con¬ 
clusion  that  Uir)  is  monotonically  increasing. 

APPENDIX  B:  ASYMPTOTIC  APPROXIMATION  OF  p(r) 

To  derive  an  asymptotic  approximation  for  pir)  as 
r— >  -1-00,  we  need  to  deal  with  integrals  such  as 
/o’"  exp(r  cos  29)d9.  First,  we  use  the  symmetry  of  cos  29  to  Combining  Eqs.  (B3)  and  (B5)  gives  us 


(B4) 


:4)  exp[r(l  -2j)]  ^  r— 
>0  V^Vl-^ 


=  4  exp(r) 


J  exp(-  2rs) 


1  H-  -S  -I- 
2 


ds 


^  ,,/r(3/2)  1  r(5/2) 

=  4  exp(r)  i=  +  -  ■ - — i=  + 

\  2rV^  2  (20^/2^ 

1  3  1 


(B5) 


rewrite  the  integral  as 

r27r 


'‘tt/I 


rzTT  rr 

I  exp(r  COS  26^)^/^=  4  I  exp(r cos  (Bl) 

Jo  Jo 

38 

To  apply  Watson’s  lemma  on  the  integral,  we  use  a  change 
of  variable, 

5  =  (1  -  COS  29)12, 

ds  =  sin  29d9=  2  sin  6cos  9d9=2^s^!l  -  sd9,  (B2) 


1  -  (cos  29)  = 


/o’’(l  -  cos  20exp(r  cos  29)d9 


/o’"  exp(r  cos  29)d9 
113  1 

2  r  16  1  /  1 

- =  —  l-F  — H- 

1  1  2r\  4r 

1  +  -.-  +  ... 

8  r 


(B6) 


and 


1  1 
d9=  — 


=ds. 


2.  •is'll' 

Applying  the  change  of  variables  to  the  integral  yields 


(cos  29)  =  1  -  (1  -  (cos  20)  =  1  -  —  H- 

2r 


(B7) 


To  find  an  asymptotic  approximation  for  function  pir),  we 
rewrite  it  as 
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p{r): 


1 


(i-) 


1+1 


(cos20)\g(r) 

1  1  1  -(cos2^) 

{cos20)  g(r)  (cos  26*) 


(B8) 


The  first  term  on  the  right-hand  side  is  exponentially  small. 
The  second  term  can  be  expanded  using  results  (B6)  and 
(B7)  obtained  above, 


p(r)=- 


1  -  (cos  20} 
(cos  20} 


1  /  1 
—  1  +  —  + 
2r\  4r 


+  ■■■ 


-if  3 
I  1  H-  H-  •  •  • 
2r  \  4r 


(B9) 


Now  we  derive  the  expansion  of  function  p(r)  in  terms  of  [/. 
As  +°°,  U,  as  a  function  of  r,  has  the  expansion 


U{r)  = 


r 

(cos  20} 


1 

l-F  — 
2r 


+  ■■■  . 


(BIO) 


Expressing  r  in  terms  of  U  yields 


r{U)  =  U\ 


+  ■■■  . 


(BID 


Substituting  into  Eq.  (B7),  we  arrive  at 


p{U)  =  p{r{U))  ■■ 


-1  / 


iH- 


2riU)\  4riU) 


-I  5 

—  1  +  —  H- 

2U\  4U 


(B12) 
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